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§1. Introduction. 


This report presents an analysis of the structure of transient 
acoustic waves in stratified fluids whose densities and sound speeds are 
functions of a single depth coordinate. The acoustic field in such a 
fluid may be described by a real-valued function u(t,x,y) (the acoustic 
potential or the excess pressure) that satisfies the wave equation 

(1.1) |^t - c 2 (y) p(y) V • (p'VyJVu) - f(t,x,y) 

where t is a time coordinate, y is a depth coordinate, x ■ (x lf x 2 ) are 
Cartesian coordinates in a horizontal plane, V * (3/3x t ,3/3x 2 ,3/3y) and 
f(t,x,y) is a function that characterizes the wave sources. c(y) and 
p(y) are the variable sound speed and density, respectively, and 
p' 1 (y) - l/p(y). 

This report is a sequel to the author's report [12] on "Spectral 
Analysis of Sound Propagation in Stratified Fluids." That work contains 
a spectral analysis of the acoustic propagator 

(1.2) Au - -c z (y) p(y) V • (p ‘(y)7u) 

for the cases of an unlimited fluid (-« < Xj,x 2 ,y <<*>), a semi-infinite 
layer (-« < x 1 ,x 2 < « and 0 < y < “) and a finite layer (-» < Xj,x 2 < <*> 
and 0 < y < h < <*>). The integration of (1.1) below is based on the 
spectral -analysis of [12] and the notation and results of [12] are used 
throughout this report. As in [12], only the case of an unlimited fluid 
is presented in detail. The modifications required in the second and 
third cases are described in §6 at the end of the report. 


1 





» 




t 


t 


f 






2 


The analysis of the structure of transient acoustic waves in 
plane stratified media was initiated by the author in 1973 [7]. This 
work dealt with the special case of the Pekeris profile (propagation in 
a half-space, p(y) - const., c(y) piecewise constant). In 1974 the 
results were extended to the symmetric Epstein profile (p(y) * const., 
c 2 (y) ■ K sech 2 (y/H) + M) [8]. The general Epstein profile 
(p(y) ■ const., c “(y) * K sech 2 (y/H) + L tanh (y/H) + M) was treated 
in 1979 [11] using the spectral analysis of the Epstein operator due to 
Guillot and Wilcox [2, 3]. A preliminary version of [11] was announced 
in 1978 [1]. The analysis of these cases was based on explicit 
representations of their normal mode functions by means of well-known 
special functions. Extension of the analysis to larger classes of strat¬ 
ified fluids had to await the extension of the spectral analysis of the 
acoustic propagators (1.2) to such classes. 

The class of stratified fluids treated in [12] was characterized 
by the properties 


(1.3) 

(1.4) 


(1.5) 


p(y) and c(y) are Lebesgue measurable, 

0 < p m < P(y) < P M < °°> 0 < c m * c( y> s «M < 



p(±«)| 


dy < », 


+00 

| c(y) - c (±°°) I dy < », 

Jo 


where p^, Pj^, p(±°°), c m> c^ and c(±°°) are constants. In this report the 
results of [10] and [11] on transient acoustic waves in stratified fluids 
are extended to the class of fluids characterized by (1.3), (1.4) and 
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( 1 . 6 ) 


|p(y) - P(±°°) | < C(±y) 


. 

|c(y) - c(±<=°)1 < CCiy) -0 

j 


for ±y > 0, 


where C and a are constants and 


(1.7) a > | . 

It is clear that (1.6) implies (1.5). It will be seen from the analysis 
below that (1.6) could be replaced by other order conditions at y * ±°°. 

It is not known whether the results of this report hold for the entire 
class of luids defined by (1.3), (1.4), (1.5). 

The remainder of this report is organized as follows. In §2 the 
normal mode expansions for unlimited fluids of [12] are used to decompose 
transient acoustic fields into free and guided components and to obtain 
integral representations of these components. The behavior for large t 
of the free and guided components is derived in §3 and §4, respectively. 
In §5 these results are used to calculate the asymptotic distribution for 
large t of the wave energy. §6 presents corresponding results for 
semi-infinite and finite layers. 








§2. Normal Mode Expansions of Transient Acoustic Fields . 

The integration of the wave equation (1.1) will be based on the 
spectral analysis of the acoustic propagator A in the Hilbert space 
3C - L 2 (R 3 ,c 2 (y) p *(y) dxdy), as developed in [12]. To this end (1.1) 
is interpreted as the equation 

(2.1) 4tt + Au = f(t,*)» t e R, 


for a function t *► u(t,*) € 3C. The wave sources will be assumed to act 
during a time interval [0,T], so that supp f C [0,T]. The corresponding 
acoustic wave is the solution of (2.1) that satisfies the initial 
condition 

(2.2) u(t,») ■ 0 for all t < 0. 

The solution, based on the spectral theorem for the selfadjoint operator 
A > 0, is given by Duhamel's integral 

(2.3) u(t,*) - | {A 1/2 sin (t - x)A l/2 } f(T,*)dt, t > 0. 

* o 

Indeed, if f € C([0,T],JC) then (2.3) is the unique "solution with finite 
energy" of [6], while if f 6 C([0,T],D(A 1/2 )) then (2.3) is the "strict 
solution with finite energy." In addition, if f € C([0,T],D(A _1/Z )) then 

(2.4) u(t,x,y) * Re (v(t,x,y)} 
where v(t,*) is the complex-valued potential defined by 
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(2.5) v(t, •) - i exp {-itA l/2 } A~ 1/2 j exp {ixA 1/2 } f(t,*)dt. 


In particular. 


( 2 . 6 ) 


v(t, •) » exp {-it A 1 ^ 2 } h for all t > T 


where 


(2.7) 


X 

i A ^ 2 | exp {it A 1 / 2 } f(t,*) dt. 


The initial value problem 


( 2 . 8 ) 



+ Au 


0 for t > 0, 


(2.9) 


u(0) 


f du(0) 
’ dt 


g 


can be treated by the same formalism. Indeed, if f 6 D(A ly<2 ) and 
g e D(A -1 ^ 2 ) then the solution of (2.8), (2.9) is given by (2.4), (2.6) 
with h * f + i A _l/2 g € D(A l/2 ); cf. {9, Ch. 3], 

The integral 

(2.10) E(u,K, t) - | { )Vu| 2 + c” 2 (y)||^| } p -1 (y) dxdy 

may be interpreted as the energy of the acoustic field u in the set 
K C R 3 at time t. Moreover, A is the selfadjoint operator in K asso¬ 
ciated with the sesquilinear form A on JC defined by D(A) - Lj(R 3 ) C Jf 


and 
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( 2 . 11 ) 


A(u,v) 


Vu • 7v p 1 (y) dxdy 

Jr 3 


It follows from Kato's second representation theorem [4, p. 331] that 
D(A 1/2 ) * L 2 (R 3 ) and for all u € D(A l/z ) 


(2.12) I)A 1 ^ 2 uflj^ * A(u,u) = |Vu| 2 p 1 (y) dxdy. 

Jr 3 

Hence the total energy satisfies 


(2.13) 


E(u,R 3 ,t) = DA 172 ulj 


3u 

at 



If h € D(A l7z ) and u is defined by (2.4), (2.6) then a simple calculation 
shows that 

(2.14) E(u,R 3 ,t) ■ BA 1/2 hB^. for all t > T. 

The analysis of the structure of the acoustic potential (2.4), 
(2.6) presented below is based on the normal mode expansion for A of 
[12, §8]. The orthogonal projections {P + ,P_,P 0 ,P 1 ,P 2 ,•••} in 3C defined 
by the normal modes form a complete family that reduces A [12, 
Corollaries 8.4 and 8.10]. Hence the same is true of the orthogonal 
projections 


(2.15) 


P f =* P + + P_ + P 0 


and 


P 

g 


No-1 

l V 

k-1 


(2.16) 




This provides a decomposition 


(2.17) u(t,*) = u f (t,*) + u g (t,») 
into orthogonal partial waves 

(2.18) u f (t,‘) = P f u(t,-)» u g (t,*) = P g u(t,*). 

The function u^, called the free component of u, will be shown to behave 
for large times like a diverging spherical wave. The function u , called 

o 

the guided component of u, will be shown to be approximately localized 
near one or more of the planes y * const, where c(y) has local minima and 
to propagate outward in horizontal planes like a diverging cylindrical 
wave. This second component shows the profound effect of acoustic ducts 
on transient acoustic waves. It is absent when c(y) has no local 


minima.. 



§3. Transient Free Waves. 


The normal mode expansions of [12] are used in this section to 
calculate the asymptotic behavior for t ■+• “ of the free component 
u^(t,*) * u(t,•) • The principal result is that in each of the half¬ 
spaces R_|_(d) and R*(d), where R 3 (d) - {(x,y) : ±(y - d) >0}, u f (t,*) is 
asymptotically equal to a wave function for a homogeneous fluid with 
propagation speed c(°°) and ci- 00 ), respectively. It is this behavior 
that motivates the term "free component" for u^t,*). 

It will be assumed that the total acoustic potential u satisfies 

u(t, *) * Re {v(t, •)} where v(t,*) ■ exp {-it A^ z } h and he D(A 1/2 ) 

(see §2). The corresponding partial waves u^(t,*) * u(t,*) with k > 1 

satisfy u k (t,*) * Re (exp (-it A 1/2 ) P^h}. This follows from the fact 

that the normal mode functions (y,p) are real for k > 1, which implies 

that P k (h) * P k (h) . It follows by addition that u^(t, •) * 

Re {exp (-it A 1 / 2 ) P h} * Re {v (t,*)} and hence 
8 8 

[ u f (t,*) - Re {v f (t,*)}. 


(3.1) 


v f (t,*) - exp (-it A l/2 ) P f h - exp (-it A 1/2 ) h f . 


The normal mode representations 


(3.2) v f (t,x,y) - <j> + (x,y,p,q) exp {-itl 1/2 (p,q)} h + (p,q) dpdq 

JR 3 ~ 

of [12, §8] provide the starting point for calculating the asymptotic 
behavior of u^t.x,y) for large t. The integral in (3.2) will, in 
general, converge only in X. For brevity the JC-lim notation of [12] is 
suppressed in this report. 
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Equation (3.2) gives two representations of v^ corresponding to 
the two families <J> + and <p_. The calculations below are based on the 
^-representation which has been found to yield the simplest form of 
the asymptotic wave function. It will be convenient to introduce the 
characteristic functions x+» Xo and x_ of the cones C + , C 0 and C_ in 
(p»q)~space [12, (1.48)—(1.51)] and to decompose h_ as 

(3.3) h_(p,q) = £(p,q) + m(p,q) + n(p,q) 

/V /V A 

where Z ■ x+ h_, m * Xo h_ and n ■ x_ h_. The corresponding decomposi¬ 
tion of v^ is 

(3.4) v f “ v £ + v m + v n 
where 

* 

v^ - exp (-it A 1/>2 ) $>*& 

(3.5) ■ v ■ exp (-it A 1 ^ 2 ) $*m 

m — 

v^ * exp (-it A 1 ^ 2 ) $*n 

The behavior for t *► 00 of these three functions will be analyzed 
separately. 

Behavior of v^. The partial wave v^ has the representation 

(3.6) v Jl (t,x,y) - J c 4>.(x»y»P»q) exp(-it w + (p,q)) £(p,q) dpdq 
where 


(3.7) 









(Recall X(p,q) * w^(p,q) for ±q > 0.) To discover the behavior of 
v^(t,x,y) for (x,y) G R 2 (d) and t -*■ 00 it will be convenient to write 
<J>_(x,y,p,q) in a way that puts in evidence its behavior for y -*• ±°°. To 
this end recall that by [12, (1.53)-(l.60)] 

(3.8) <J>_(x,y,p,q) = (2 tt)“ 1 c(») (2q) 1/2 e ip X ^ + (y» |p [» ^) 

for (p,q) = X + (p,X) € C + (and hence X = X(p,q) ■ c 2 (“)(|p| 2 + q 2 )). 
Moreover, by [12, (4.5), (4.6) and (8.41)] one can write 


(3.9) 


, 1/2 

( CO ] 1/2 

“ l4<(l-)j T + (y,V,\) exp {-iy q_(y,X)} 


where 


(3,10) T + (y,y,X) - T + (y,X) 4\(y,U,X) exp {iy q_(y,X)}-*-T + (y,X) , y -*■ 

Similarly, by [12, (4.1) and (4.6)-(4.11)] one can write 


(3.11) 


, 1/2 

4> + (y,V,\) - I I + <y»M.X) exp {-iy q + (u,X)} 

+ R + (y»P»X) exp (iyq + (y,X)}] 


where 

' ^ 
i + (y»u,X) - 4> 2 (y,y,X) exp Uyq + (y,X)} -*■ l 

(3.12)' • y -*■ +». 

R + (y»y»X) - R + (U,X) <h(y,u,X) exp {-iy q + (y,X)}->- R + (y,X) 


Combining (3.6), (3.8) and (3.11) gives 
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VjU.x.y) 


c(°°) p(°°) 1/2 (2 tt)~ 3/1,2 


exp {i(x • p+yq - tw + (p,q) } I + (y, |p j, X) £.(p,q)dpdq 


(3.13) 


+ c(°°) P(°°) 1/2 (21T) 3/2 j exp (i(x * p-yq - tu + (p,q) } R + (y, |p |, X) £(p,q)dpdq. 

* r 


It is natural to expect that in R^(d) the partial wave v^(t,x,y) will 
propagate as t + ® into regions where y is large and hence I + (y, |p |, X) 
and R_j_(y, Ip |»X) are near their limiting values. Thus the representation 

(3.13) suggests the conjecture that 

(3.14) v^t,*) ~ +v^(t,«) inL 2 (R 3 (d)), t ■+ °° 

where v£ and v^ are defined by 


(3.15) v £(t,x,y) = c(°°)p(°°) (2 tt) 3 ^ 2 j exp {i(x*p+y q - tu> + (p,q) K(p,q)dpdq 


v£(t,x,y) - c(°°)p(®) 1/2 (2 tt)“ 3/2 j exp (i(x*p-yq - tu> + (p,q)} x 


(3.16) 


x VlPl> A) *(P»Odpdq 

- c(»)p(») 1/2 (27T)~ 3/2 | exp {i(x*p+y q - tui + (p,q)} x 

-C + 

x R + (|p|,A) £(p,-q)dpdq 


where -C + - {(p,q) : (p,-q) € C + } - {(p,q) : q < -a|p|}. Note that v® 
and v^ are waves in a homogeneous medium with density p (<*>) and sound 
speed c(«°) . More precisely, 
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(3.17) 


vUt, *) 


v*(t,-) 


exp (-i t c(“) Aj /2 ) h c 


exp (-i t c(“) k 1,z ) hi 


where A 0 is the selfadjoint realization in L 2 (R 3 ) of 

-A ■ -(3 2 /3xJ + 3 2 /3x 2 + 3 2 /3y 2 ) and h^ and h^ are the functions in 

L 2 (R 3 ) whose Fourier transforms are 


(3.18) 


h^p.q) * c(“) p(“) 1/2 i-(p,q) - c(°°) p(°°) U2 X^(p,q) h_(p,q), 


1/2 


h^(p,q) - c(«) p(°°) 1/2 R + (|pj,X) 2(p,-q) 


c (») p(») 1/2 R + (|pj,X) (1 - X+(P»q)) h_(p,-q) 


Both functions are in L 2 (R 3 ) because h_ € L 2 (R 3 ) and lR + (|p,A)| <_p 1/,2 («>) by 

A t 

the conservation law [12, (1.46)]. Moreover, supp h^ C -C + and hence 
the theory of asymptotic wave functions for d'Alembert's equation 
[9, Ch. 2] implies that v^(t, •) - 0 in L 2 (R+(d)) when t -*• 00 . Combining 
this with (3.14) gives 


(3.19) v Jl (t,*) ~ v£(t,*) in L 2 (R|(d)), t ■+ °°. 

Now consider the behavior of v^(t,x,y) for (x,y) e R 3 (d), t -*■ °°. 
Combining (3.6), (3.9) and (3.10) gives 


VjU,x,y) - c(«) p(“) 1/2 (2 tt) -3/2 j exp {i(x*p + yq_ - tu> + ) } x 


(3.20) 


* T + (y, |p ,A) £.(p,q)dpdq 


where q_ - q_(|p|,A), w + - w + (p,q) and A - A(p,q) - u>+(p,q). This 
representation suggests that 
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(3.21) 


v 5 (t,*) ~ v»(t,*) in L 2 (Rj(d)), t -*• “ 


where 


(3.22) 


v|(t,x,y) - c(») p(°°) 1/2 (2 tt)“ 3/2 j exp (i(x • p + yq_- tco + )} * 


r + ( i P | »*) *-(p,q)dpdq 


Now the mapping (p,q) ■+• (p.q 1 ) =* X'(p,q) ■ (p»q_( |p| »^(p,q))) with 
domain C + has range X'(C + ) = Rj[ ■ R+(0), Jacobian 9(p,q)/3(p,q') 
m c 2 (-°°)q , /c 2 ( 00 )q and satisfies w + (p>q) = oj_(p,q’). Thus (3.22) implies 


the representation 


Vo(t,x,y) - c(°°) p(“) l/2 (2 tt) -3/2 f exp (i(x ■ p + yq' -tu> )} * 

* Jr* 


(3.23) 


* T + (|p| , 0 J^(p,q')) 2-(p, q) (c 2 (- 00 ) q '/c 2 (°°)q) dpdq ' 


where q ■ q(|p|,q') * /a 2 ( p 2 + q' 2 ) + q' 2 . Note that 


(3.24) 


v 2 (t,-) - exp (-itc(-») Aj /2 ) h| 


where h| e L 2 (R 3 ) has Fourier transform 


(3.25) h^(p.q') - c(«) p(») 1/2 T + ( |p| ,U 2 (p,q')) Jl(p,q) (c 2 (-«)q 7c 2 («)q) . 

Since supp h£ C the results of [9, Ch. 2] imply that v|(t,*) ~ 0 in 
L 2 (R 3 (d)) when t ■+• 00 . Combining this with (3.21) gives 


(3.26) 


v 9 (t,*) - 0 in L 2 (R 3 (d)), t -*• «. 


Analogous conjectures concerning v m (t,*) and v n (t,*) will now be 
formulated. Only the main steps of the calculations will be given since 
the method is the same as for v«(t,»). 
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Behavior of v . v has Che representation 
_ m m 

(3.27) v m (t,x,y) = 4>_(x,y,p,q) exp (-1 1 U) + (p,q)) m(p,q)dpdq, 

Co 

by (3.5), where 

(3.28) <j>_(x,y,p,q) - (2ir) -l c(°°) (2q) e ip ‘ x i|» 0 (y, |p| ,X) 

for (p,q) ■ X 0 (p,X) 6 C 0 (and hence X « X(p,q) ■ u>£(p,q)). Moreover 
(see [ 12 , <4.18)-(4.24)]), 


cfoo} } 1/2 

( 3 . 29 ) 'P 0 (y,u,x) - [4 ~ ffqT(y, ~ i 7 j T o(y»i*» x > ex P (yql(vi>*>) 


where 


(3.30) T 0 (y,y,X) - T 0 (y,X) <j> 3 (y,y,X) exp (-yq^(y.X)) -*• T 0 (y,X), y -*■ -» 


and 


(3.31) 


where 


^o ( y»M> “ ^iffiy.i) ] ' UoCy.M.A) exp {-iy q + (y,X)} 
+ R 0 (y.y»*) exp {iyq + (y,X))J 


(3.32) \ 


I„ (y»y*X) - 4> 2 (y,y»X) exp {iy q + (y,X)} -*■ 1 


R 0 (y.M.X) - R„(y,X) <()j(y,y,X) exp {-iy q + (y,X)}-^ R 0 (y,X) 
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Combining (3.27), (3.28) and (3.29) gives 


v m (t,x,y) = c(°°) pO 30 ) 1 '* 2 (2 tt) 3/2 


exp (i(x • p - tw + )} x 


(3.33) 


xT 0 (y,|p|,X) exp (y q^) m(p,q)dpdq 


where u) + * oo + (p,q), X * uj^(p,q) and q^ * q_^C |p|»X>. Since 


(3.34) 


T (y, |pI.X) exp (y q_l(|p|.X)) -*-0, y 


equation (3.33) suggests that 


(3.35) 


v m (t,*) - 0 in L 2 (R 3 (d)), t -*■ ®. 


Similarly, combining (3.27), (3.28) and (3.31) gives 


v m (t,x,y) 


- c(°°) p(«°) 1/z (2ir) 


»/2 (">^" 3/2 


exp {i(x • p + y q - tw + )} x 


x I o<y>|p|* x > m (p,q)dpdq 


(3.36) 


+ c(°°) p(°°) 1/2 (27r) 3/2 | exp {i(x • p-yq - tu) + )> 


x R„(y,Ip ,X) m(p,q)dpdq 


which suggests that 


(3.37) v m (t,-) ~ v®(t,*) + v 3 (t,*) in L 2 (R 3 (d)), t - ® 


where v® and v 3 are defined by 
mm 


(3.38) 


v°(t,») ■ exp (-itc(») A 3/2 ) h ffl 


v 3 (t,*) - exp (-itc(») Ai /2 ) h 3 


> 
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and h and h l are the functions in L,(R 3 ) whose Fourier transforms are 
mm 

h m (p.q) * c(«) p(®) 1/2 m(p,q) - c(») p(«>) 1/2 Xo(P»q) h_(p,q), 

(3.39) ■ 

h^(p.q) * c(») p(oo) 1 /2 r o (| p| f A) m(p,-q) 

- c(°°) p(“) 1/2 R 0 (|p| ,X) (1 - X 0 (p»q)) h_(p,-q) . 

Note that supp h 1 C R 3 and hence v 1 (t,*) - 0 in L,(R?(d)) when t •*> 

* m - in *• “r 

Combining this with (3.37) gives 

(3.40) v m (t, *) " v£(t, •) in L 2 (R^(d)), t -*• <». 

Behavior of v . v has the representation 
_ n n 

(3.41) v (t,x,y) * ^(x,y,p,q) exp (-itoj (p,q)) n(p,q)dpdq, 

'C_ 

by (3.5), where 

(3.42) <j>_(x,y,p,q) - (2TT) _ 1 c(-«)(2|q |) 1/2 e ip ‘ x 4<_(y, | P | ,X) 

for (p,q) ■ X_(p,X) € C_ (and hence X ■ X(p,q) * to 2 (p,q). Moreover 
(see [12, (4.5) and (4.12)-(4.17)]), 

(3.43) ^_(y,U»X) - [ 4 ^^X) ] T_(y,p,X) exp Uyq + (y,X)} 

where 

(3.44) T_(y,y,X) - T_(y,X) <J>, (y,y,X) exp {-iy q + (y,X)} -* T_(u,X), y -*• « 
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and 


(3.45) 


i|;_(y,p,X) 


. P(-°°) 1 1/2 

4irq_(p,X)J 


[I_(y,P»X) exp (iy q_(p,X)} 


+ R_(y,U,X) exp {-iy q_((i,X)}] 


where 


I_(y.U,A) - 4> 3 (y,U,X) exp {-iyq_(p,X)} -*• 1 


(3.46) < 



R _(y.V»*) - R_(V,X) $ 4 (y»V»X) exp {iy q_(y,X) } -*• R_(v* X) 


Combining (3.41), (3.42) and (3.45) gives, after simplification using 
q_(|p|» ( * J I(p»q)) - (q 2 ) 12 ° -q for (p,q) € c_. 


v n (t,x,y) 


(3.47) 


c(-«) p(-“) (21T) 3 ^ 2 exp {i(x • p+yq - tin )} * 

J c 


* x _(y.|P .*) n(p,q)dpdq 


+ c(-») p (-°°) 1,1 (2 tt)~ 3/2 exp{i(x* p-yq - to) )} 

■*C_ 

x R _(y»|P1»*) n(p,q)dpdq. 


This suggests the asymptotic behavior 

(3.48) v Q (t,*) ~ v^(t, •) + v^(t,«) in L 2 (R 3 (d)), t ->• ®, 


where 


(3.49) 


v®(t,*) - exp (-itc(-®) A^ /2 ) h Q 
v^(t,*) - exp (-i t c(-°°) A^ 2 ) h^ 
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and h and h 1 are the functions whose Fourier transforms are 
n n 


(3.50) 


h n ( p ,q) = c(-°°) p(-°°) 1/2 n(p,q) = c(-“) p(-°°) 1/2 X_(P>9) h_(p,q). 


| h^(p,q) * c(-°°) p(-°°) 1/2 R_(|p| ,X) n(p,-q) 

= c(“°°) p(-oo) !/ 2 R_(|p|,X) (1 - X_(P.9>) h_(p»-q). 

Note that supp h^ C and hence v^(t,*) ~ 0 in L 2 (R^(d)) when t ■+• 00 . 
Combining this with (3.48) gives 

(3.51) v n (t ** ) - v n (t ’ #) i* 1 L 2 (R 2 (d)), t ♦ ». 

Finally, combining (3.41), (3.42) and (3.43) gives 


v n (t,x,y) 


= c(-°°) p(-°°) 1/2 (2tt)~ 3/2 


exp {i(x • p - y q + - tai_) } x 


(3.52) 


x T_(y,Jp| ,X> n(p,q)dpdq 

where w_ * w_(p,q), X =» w 2 (p,q) and q_j_ ■ q + (|pj*X). This suggests that 
(3.53) v n (t >*) ~ v n (c ** ) i*»L 2 (R*(d)), t » 


where 


v 2 (t,x,y) 

(3.54) 


c(-oo) p (-«) 1/2 ( 27T) “ 3/2 


exp (i(x • p - y q + - tu>_ 


)} 


x T _(|p|,X) n(p,q)dpdq. 


Now the mapping (p,q) (p,q') - X"(p,q) - (p,-q + (p,X(p,q))) maps C_ 

onto X"(C_) ■ -C + , has Jacobian 3(p,q)/3(p,q') ■ c 2 (°°)q'/c 2 (-°°)q and 
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satisfies w_(p,q) *u> + (p,q'). Thus 


(3.55) 


v*(t,‘) - exp (-i t c(°°) Aj /2 ) h 2 


where h 2 has the Fourier transform 
n 


h*(p,q') * c(-“) P(-°°) 1/2 T_(|p[,w£( p ,q')) n(p,q(p,q')) 


(3.56) 


x (c 2 (°°)q '/c 2 (-°°)q) 


Moreover, supp h* C R 3 and hence v^(t,*) ~ 0 in L 2 (R+(d)), t 
Combining this with (3.53) gives 


(3.57) 


v n (t.*) ~ 0 in L 2 (R^(d)), t -*• “. 


The asymptotic behavior of v^(t,*) for t -*• 00 may be obtained from 
the three cases analyzed above by superposition, equation (3.4). Thus 
equations (3.19), (3.26), (3.35), (3.40), (3.51) and (3.57) imply 


v£(t,•) + v£(t,*) in L z (R^(d)) 


(3.58) 


v f (t,*) 


v“(t,*) 


in L,(Rf(d)) 


t <*>. 


On combining this with the definitions of v£, v£ and v£, equations (3.17), 
(3.18), (3.38), (3.39), (3.49) and (3.50), one is led to formulate 
Theorem 3.1 . For every h 6 K let v£(t,») be defined by 


(3.59) v°(t,x,y) 


[exp (-itc(«) Aj /2 ) h + (x,y), (x,y) € R^.(d), 


exp (-itcH A, /2 ) h“(x,y), (x,y) € Rf (d), 


where h and h~ are the functions in L 2 (R 3 ) whose Fourier transforms are 


given by 




i 
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(3.60) 


h + (p,q) 


c (oo) p(«) 1/2 h_(p,q) , (p,q) € R+, 


, (p,q) € R 3 , 


and 


(3.61) 


_ 

h (p,q) 


0 , (p,q) £ R+, 

c (_o°) p(-°°) 1/2 h_(p,q) , (p,q) SR 3 . 


Then 


(3.62) 


lim 8v f (t,•) - v° (t,*)8 

fc-XJO 


X 


0 . 


Theorem 3.1 implies corresponding asymptotic estimates for the 
free component u f (t,*) = P f u(t,‘) = Re {v f (t,*)} of the acoustic potential 
u(t,*). Indeed, if u^(t,*) is defined by 


(3.63) u°(t,*) » Re (v£(t,‘)} 

then Theorem 3.1 and the elementary inequality |Re z| < |zj imply 
Corollary 3.2 . For all h € X one has 


(3.64) 


lim 8 u,(t, •) 

t-KO 


Uf (t, = 0. 


If the initial state h has derivatives in Jf then u^(t,*) and 
u°(t,*) have the same derivatives in X and (3.64) can be strengthened to 
include these derivatives. In particular t one has 

Corollary 3.3 . For all h € L*(R 3 ) - D(A 1/2 ) one has 

(3.65) lim ||D, u f (t,*) - D. u£(t,*)8jf - 0, j = 0,1,2,3, 


J 
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where D 0 * 3/3t, Dj “ 3/3x lt D 2 * 3/3x 2 and D 3 = 3/3y. Equation (3.65) 
is equivalent to convergence in energy: 

(3.66) lim E(u^ - u£,R 3 ,t) = 0. 

t-x*> 

Corollary 3.3 can be proved by applying the method of this 
section to the derivatives v^(t,*) (j ■ 0,1,2,3) which are given by 

integrals of the same form as (3.2). Detailed proofs for the case of the 
Pekeris profile were given in (10]. 

Proof of Theorem 3.1 . The remainder of this section is devoted 
to the proof of Theorem 3.1. The decomposition (3.3) is used for the 
proof. Moreover, for brevity, only the asymptotic equality (3.53) for 
v n (t,*) is proved. The remaining five cases, namely (3.14), (3.21), 
(3.35), (3.40) and (3.48) can be proved by the method used for (3.53). 

As a first step, (3.53) will be proved for the special case of n(p,q) 

€ C 0 (CJ, the set of continuous function with compact supports in the 
open cone C_. The general case will then be proved by using the fact 
that C 0 (C_) is dense in L 2 (C_). 

For functions n(p,q) € C 9 (C_) the integrals defining v^ and v* 
converge point-wise, as well as in JC, and one can write 



where 





23 


(3.68) w(y,p,t) 


ro 


exp (-i(yq + + tu>_)} [T_(y, |p|,o>*) - T_( jp|,cof) ]n(p,q)dq 


Parseval's formula in L 2 (R 2 ), applied to (3.67), gives 


(3.69) [ |v (t,x,y) - v 2 (t,x,y) | 2 dx«c 2 (-«) p(-“») (2tr)~ 1 ( |w(y,p, t) | 2 dp. 

r 2 n n Jr 2 


On integrating this over y > d one finds 


(3.70) lv n (t,*) - v 2 (t,‘)l L2(R 3 (d)) -c(-») p(-°°) 1/2 (2 tt )~ 1 ^ 2 lw(*,t)l T /d3 


L 2 (R^(d)) 


The last relation implies that to prove (3.53) it is sufficient 
to prove that 


(3.71) 


w(*,t) + 0 in L 2 (R|(d)), t -*■ “>. 


To this end it will be convenient to change the variable of integration 
in (3.68) from q to u ■ oa_(p,q) ■ c(-°°) /|p | 2 + q 2 . Solving this 
equation for q < 0 gives q ■ -(d) 2 c 2 (-°°) - |p| 2 )^ 2 * ~<5_( |p| »u) 2 ) with 
oj > c(-°°)|pj. Hence (3.68) can be written 


(3.72) 


w(y,p,t) = exp (-.11 u>) W(y,p,u>)du) 

) c (-®)|pj 


where 


W(y,p,ai) 

(3.73) 

c 2 (-°°) exp {-iy q^.( |p | ,o> 2 ) }[T_(y, |p | .u)*) - T_( jp | ,cj 2 ) ]n(p,-q_( |p| ,ou 2 ))u) 

q_(|pj ,co 2 ) 
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The assumption that n e C 0 (C_), together with (3.44), implies that 
W £ C 0 (R x O where T » {(p,w) : <u > c(-°°)|p|}. Moreover, by a standard 
partition of unity argument, one may assume without loss of generality 
that 

(3.74) supp W(y,*) C {(p,w) : |p| < p 0 and 0 < u) 0 < w < Wj} 
for all y € R where oj o > c(-«)p a . This in turn implies that 

(3.75) w(y,p,t) - I 1 exp (-itu) W(y,p,u)daj 

H 

and 

(3.76) supp w(y,•,t) C B(p 0 ) - (p : |p| < p Q } 
for all y € R and t € R. Thus 


lv(«,t)|f ( j,.., - I [ |w(y,p,t)( 2 dpdy 


(3.77) 


[ [ |w(y»P,t )| 2 dpdy 

J d 'B(p 0 ) 


+ j J |w(y,p,t)| 2 dpdy 

J y 0 ^( Po ) 


for any y fl > d. The proof of (3.71) will be derived from (3.77) and the 
following two lemmas. 

Lemma 3.4 . Let n€ C„(C_) and assume that (3.74) holds. Then 
for each d € R, y fl > d and p 0 > 0 one has 
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(3.78) lim w(y,p,t) - 0, 

t-*» 

uniformly for all (y,p) S [d,y 0 ] x B(p 0 ). 

Lemma 3.5 . Under the hypotheses of Lemma 3.4 there is a constant 
C * C(n) such that 

(3.79) |w(y,p,t)| < C y 1_a 

for all y > 0, p £ B(p 0 ) and t e R, where a > 3/2 is the constant of 
condition (1.6). 

Proof of Lemma 3.4 . The proof is based on a well-known proof of 
the Riemann-Lebesgue lemma. Note that by (3.75) one has 


w(y,p,t) 


(3.80) 


■<i) 1 -(7r/1) 

exp (-iw t) W(y,p,art-(7r/t))du> 

^u 0 -(ir/t) 



exp (-i t w) [W(y,p,co) - W(y,p,w+(iT/t)) ]doj 


1 f^ 0 

-K exp (-itaj) W(y,p,uH-(iT/t))du) 

■’oi 0 -(tt/ t) 


+ 4 ] exp (-i t u) W(y,p,w+(iT/t))dw. 

•’cj 1 -(Tr/t) 

The limit relation (3.78) is obvious from (3.80) and the continuity of 
W. The uniformity of the limit follows from (3.80) and the uniform 
continuity of W on compact subsets of R * T. 

Proof of Lemma 3.5 . Note that by (3.72), (3.73) and (3.44) one 


has the estimate 






I 
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(3.81) 


|w(y,p,t) 



(y. |p|.<*> 2 ) exp (-iy q+(|p|»w 2 )}- l| M(p,U))dw 


for all y, p and t where 


(3.82) 


m(p,oj) - c _a (-“) |t_(| p|, u> 2 ) | 


n(p,-q_( |p|,0} 2 ))(J 


q_( IpI.w 2 ) 


is continuous for |p| < p 0 , u 0 < id < Uj. It follows that for any y € R, 

P e B(p 0 ) and t e R 

(3.83) |w(y,p,t)| < M„ sup |<J>i (y, |p| ,u> 2 ) exp {-iy q+( |p | ,w 2 ) } - l| 

where M 0 ■ M„(n) * sup M(p,u) and the suprema are taken over all |p( < p 0 
a«.d o) 0 < tii < ujj. The proof of (3.79) will be based on (3.83), conditions 
(1.3), (1.4), (1.6) on p(y) and c(y) and the proof in [12] of Theorem 
2.1; see [12, p. 27ff]. Note that because of the continuity of w it will 
suffice to prove (3.79) for all y > y t where - yi(n) is a positive 
constant. 

The solution ^(y.y.A) with X > c 2 (-®)y 2 > c 2 (°°)y 2 satisfies 
[12, p. 27ff] 

(3.84) <My,y,A) - exp {iyq + (y,A)> (n x + n 2 ) 

where n - (n 1 ,n 2 ) is characterized on y > y x as the unique solution of 
the integral equation [12, (2.64)] which can be written 


(3.85) 


n - n® + K(y,X)n, n® - (i,o). 


The kernel K(y,y',y,A) - (K^(y,y*,y,X)) is defined by 


» 
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(3.86) 


and 


K 1 : j(y,y',y,A) 


y i _ y « y 


-E 1 j(y , ,Vi,X), y' > y. 


(3.87) K (y.y'.y.A) = 


o » y x $ y' < y. 

|-exp {-2i(y-y')q + (u,X) E 2 j(y',u,X), y' > y. 


where [12, p. 27] 


(3.88) 


E(y,y,A) - B -1 (y,A) N(y,y,A) B(y,A). 


From these relations one has [12, (2.64)] 


<J>i(y,y,X) exp {~iy q + (p,X)} - 1 » Hj + n 2 - 1 


(3.89) 


0<X> *co 

-J E i j(y , ,y,X) n j (y')dy f - exp {-2i(y-y’)q + (y,A)} E 2 j(y'(y')dy' 


and hence 


2 fOO 

(3.90) |<j) 1 (y,y,X) exp {-iy q^y.A)}-l| < f j |E. k (y',y,X) | |n k (y') |dy'. 

j,k-l } y 

Using [12, (2.66)] and the continuity of B(y,A) on X > c 2 (-»)y 2 it can 
be shown that 

(3.91) lK(y,A)l < 1/2 for 0 < y < p 0 , u> 2 < X < w 2 
provided y t - y t (n) is large enough. Thus 
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(3.92) 


!nfl < 


- 1 - IKI 1 - IKI - 


< 2 


for 0 < y < p 0 , uij < X < wf. Combining this with (3.90) and [12, (2.66)] 
gives 


(3.93) 


,oo I 

l4>j(y.y.X) exp {-i y q, (y,X)} - l| <2 l |E.(y',y,A)j dy' 

iy j,k-l 

< C, j |p(y') - p(“)| dy' 

h 

+ C 2 j |c(y ') - c(“>) | dy' 

* XT 


for the same values of y and A where C x and C 2 depends only on p 0 , a) 0 
and It follows from (3.93) and (1.6) that 

(3.94) sup |<J),(y, |p|,(iJ 2 ) exp (-iyq + (|p|,w 2 )} - l| < y 1-Ct 

for all y > y x (n) where C 3 depends only on p„, u>„ and oj t (i.e., n) . 
Combining (3.94) with (3.81) gives (3.75). 

Proof of Theorem 3.1 (completed) . Lemma 3.5 implies that 

(3.95) f f |w(y,p,t) j 2 dpdy < (irp 2 C 2 /2a - 3) y 3 2a 

y 0 J b(po) 

where 3 - 2a < 0. Thus given any e > 0 there is a y 0 m y 0 (e,n) such that 


(3.96) 



B(P 0 ) 


|w(y,p,t)| 2 


dpdy < e 









for all t € R. Equation (3.77), together with Lemma 3.4 and the 
estimate (3.96), imply that 
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(3.97) 


lim sup Uw(*,t)li T z < e. 


L 2 (R*(d)) - 


Since e > 0 is arbitrary this implies (3.71), as required. 


The arguments given above, applied to v^, v m and v n , show that 
the conclusions of Theorem 3.1 hold for all h such that 
h_ € C 0 (C + u C 0 u C_). Moreover, this set is dense in L 2 (R 3 ) and hence 
C 0 (C + U C 0 u C_) is dense in JC f => P f X by [12, Cor. 8.12]. These 
facts can be used to extend (3.62) to all h 6 X because the mappings 
U(t) : X f *> L 2 (R 3 ) and U„(t) : X { ■* L 2 (R 3 ) defined by U(t)h f 
* exp (-It A 1 ^ 2 ) h^ and U 0 (t)h^ = v^(t,*) are uniformly bounded for all 
t e R (see [11, p. 32]). The density argument needed to extend (3.62) 
to all h € JC has been given in many places; see, for example, [9, Ch. 2] 
or [10, p. 260]. 
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§4. Transient Guided Waves . 

The asymptotic behavior for t -*• 00 of the guided component 
Ug(t,*) ■ P g u(t,*) is derived in this section. u^(t,*) is a sum in 
of mutually orthogonal partial waves u^Ct,*) ■ u(t,*) * Re (v^(t,*)}, 

1 < k < N 0 . The starting point for the analysis is the integral 
representation 


(4.1) 


v k (t,x,y) « | ^ k (x,y,p) exp (-itoo k (|p|) h k (p)dp 


where 


(4.2) 


h k ( p ) =* | 3 \( x »y» p ) Mx,y) c -2 (y) p -1 (y)dxdy 


and the integrals converge in 30 and L 2 (^ k ), respectively. The integral 
in (4.1) can be written 


(*•3) v k (t,x,y) » ~ | exp {i(x* p- t o» k (|p|))} </> k (y,p) h k (p)dp. 


This is an oscillatory integral that can be estimated by the method of 
stationary phase when 6 C™(ft k ). To apply the method define 


r » /t 2 +■ | x 1 


(4.4) 


t - r $ 0 , x x - r 5 lf x 2 - r 5 2 


€ - (C 0 ,€i.e 2 ) 6 S ! C R J 
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where S 2 denotes the unit sphere in R 3 . Then (4.3) takes the form 


(4.5) 


v fc (t,x,y) 



exp {ir 0 k (p,C)} g k (p»y) dp 


where 


(4.6) 


6 k (p.S) - ZtPt + C 2 P 2 - S 0 “k ( l p |) 
8 k (p»y) * (2ir) -1 i|> k (y,p) h^p) 


Estimates for large r of the integral in (4.5) are needed that are 
uniform for (£,y) in compact subsets of S 2 * R. Such estimates are 
provided by a version of the method of stationary phase due to 
M. Matsumura [5]. A form of Matsumura’s results applicable to (4.5) was 
presented in [11, Appendix], This result is applied below to estimating 

V 

The phase function 9 k has a point of stationary phase if and only 
if 

(4.7) U k (|pl)p/|p| - x/t 


where 


(4.8) U k (!p|) - u)^(|p|) 

is the group speed associated with the dispersion relation u - ci> k (|p|). 
It will be assumed, for brevity, that U k (y) is a monotone decreasing 
function that maps 0 k onto (c m ,c( so )). In this case (4.7) has a unique 
solution if |x|/t lies in the range of U k (|p|); that is 
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(4.9) 


c m < |x|/t < c(°°). 


and no solution otherwise. The solution is given by 


(4.10) 


P ■ Q V ( (/1) x/ jx 


where is the inverse function to U^. By calculating the Hessian 0£ 


one can show that 


(4.11) 


|det 0J|(p,C)| “ t 2 U. (jp|) |U'(|p|) j/|p| 


and sgn 9 k (p,£) ” 0. In particular, each point of stationary phase is 
non-degenerate and makes a contribution 


(4.12) v. (t,x,y,p) * 


jpj V2 exp {i(|x| |p|-ta) k (|p|))} V y, l p l)\(P) 
t{\(|pl)|u:(| P |)|} 1/2 


to the integral in (4.3), where p is given by (4.10). For |x|/t outside 
the interval (4.9) there is no point of stationary phase. Thus the 
stationary phase approximation to v k (t,x,y) is given by 


(4.13) 


v k (t,x,y) » x(l x |/t) v"(t,x,y,Q k (|x|/t) x/|x|) 


where X is the characteristic function of the interval (c m ,c(°°)) and one 


Theorem 4.1 . For all h S 5f such that h^ 6 C (J (fi k ) there exists a 
constant C * C k (h) such that 


(4.14) 


|v k (t,x,y) - v k (t,x,y)| < C/t J 


for all t > 0, x £ R 2 - (0) and y£ R. 




Theorem 4.1 can be proved by application of Theorems A.l and A.2 


of [11]. The proof that ^ k (y,p) has the required p-derivatives is 
lengthy but straightforward and will not be given here. If is not a 
smooth function then Theorems A.l and A.2 are not applicable and the 
estimate (4.14) may fail. However, the definitions (4.12) and (4.13) are 
meaningful for all h 6 X and one has 

Theorem 4.2 . For all h e JC, all t > 0 and k = 1, 2, 3,*** one has 

(4.15) v"(t,*)<S3C 

and 


(4.16) |v k (t,,)l JC " , ^lt , L a (O k ) " lP k h, 3C * 

00 

Moreover, the mapping t -*• v k (t,*) is continuous from R + to 3C and 

(4.17) lim lv k (t,*) - v“(t, •)*,<.- 0. 

The proofs of these properties are the same as those for the 
Pekeris profile, given in [10], and are not reproduced here. On 
defining 

(4.18) u£(t,x,y) - Re (v k (t,x,y)} 


one also has 

Corollary 4.3 . For all h € K and k-1, 2, 3,••*, 


(4.19) 


lim lu k (t,») - UfcU.OIjf 
t-»® 


0 . 


If h € Lj(R 3 ) then u^Xt,*) € Lj(R 3 ) and asymptotic wave functions 
for the first derivatives of u^ can be constructed. Indeed, if 
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€ C fl (^) then the first derivatives of v^ are given by 


D t v k (t ’ x ’ y) 
(4.20) 

Dj v k (t,x,y) 
D y v k (t,x,y) 


2tt 


2ir 


_ 1 _ 

2rr 



exp (i(x • p- t a) k (|p|))}(-iw k ( |p |))'(' k (y,p)h k (p)dp, 


exp (i(x • p - t u^dpl)) (ip..) ^ k (y,p) h k (p)dp (j = 


exp {i(x * p- t o^CIpI))} D y ^ k (y,p) h k (p)dp. 


1 , 2 ), 


These integrals have the same form as the integral (4.3) for v k> The 
corresponding asymptotic wave functions are defined by 


(4.21) 


v”(t,x,y,p) - (-i u^OpD) v k (t,x,y,p), 


v k j(t,x,y,p) » (ipj) v k (t,x,y,p) (j * 1,2), 


v k3 (t,x,y,p) • D y v"(t,x,y,p), and 

v kj(t,x,y) = X(!x|/t) v”(t,x,y,Q k (|x|/t) x/|x|) 

for j * 0,1,2,3. The analogue of Theorem 4.2 is 

Theorem 4.4 . For all h e L*(R 3 ), all t > 0 and k *» 1,2,3,••• one 

has 


(4.22) 


v k j(t,-) 6 L 2 (R ! ), j - 0,1.2,3, 


(4.23) l»* ^ 1 \ J < t -->^ 2 (R>,p-d*dy) ' 2 ,A ‘ /2h k'x- 
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and 


(4.24) 


11m ID 
t -»o 


V k (t **> " V kj (t * ,)B L, 


(R 3 ) 


0, j = 0,1,2,3. 


The proof of Theorem 4.4 la the same as for the special case of 
the Pekeris profile which was treated in detail in [10]. 

The preceding discussion was restricted to the special case 
where U^(y) is monotonic. If U^y) has a finite number of maxima and 
minima there are a corresponding number of points of stationary phase 
and the form of the asymptotic wave function is more complicated but 
still tractable. In the case of the Pekeris profile, treated in [10], 
there are two points of stationary phase. Cases that lead to 
infinitely many stationary points have not yet been encountered. They 
would require additional analysis. 




§5. Asymptotic Distributions of Energy for Large Times . 

The total energy of the acoustic field u(t,x,y), given by (2.14), 
is constant for t > T. The same is true of the partial waves u^, u g and 
u^, k-1, 2, 3,***. Moreover, it was shown in [12] that {P^Pj ,P 2 , • ••} 
is a complete family of orthogonal projections in X that reduces A. It 
follows that 


No-1 

(5.1) IA 1/2 hl£ - 1A 1/2 h f lj + I | A 1/2 1^1 .J. 

k-1 


which may be interpreted as an energy partition theorem. The partial 
energies 


(5.2) 


E (u f ,R 3 , t) = IA^ 2 h f l 2 , 

E(u k ,R 3 ,t) - IA* 72 h^.1 2 ., k-1, 2, 3,'• •, N 0 -l, 


can be calculated from the source function f(t,x,y) and the normal mode 
functions. The relationship between h and f is given by (2.7), which 
implies that 


(5.3) h_(p,q) « i X -1/2 (p,q) f_(X 1/2 (p,q),p,q) 


where 


(5.4) f_(w,p,q) 



exp (iwr) <P_(x,y,p,q) f (t,x,y) c" 2 (y) p _1 (y)dydxdT 


and similarly 
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(5.5) 


where 


\(p> - i \ l (|p|) ^k.(“k( IpI)*P> 


(5.6) £ k (»,p) - ll s exp (iu) t) ^(x.y.p) f(T,x, 


y) c *(y) P~ 1 (y)dydxdt. 


It follows chat the partial energies are given by 


(5.7) £(u,,R 3 ,t) - ( |f (X 1/2 (p,q),p,q)| 2 dpdq 

* Jb* 


and 


(5.8) 


E(u k ,R 3 ,t) 


| !^k (a) k ( lpl)’ p) l 2 dp * k * L » 2 > 


for every t > T. 

The theorems of §3 and §4 make it possible to calculate asymptotic 
distributions of energy in bounded and unbounded subsets of R 3 . Only the 
principal results are formulated here. The proofs are omitted since 
they are the same as those for the Pekeris profile which were given in [10]. 

The notation 


(5.9) E*(u,K) ■ lim E(u,K,t) 

t-K«> 

will be used whenever the limit exists. A first result is the transiency 
of all wavers with finite energy in stratified fluids: 

(5.10) E°°(u,K) - 0 for all compact sets K C R 3 . 
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The Free Component . The results of §3 imply that in each half¬ 
space R^(d) (resp. R 3 (d)) behaves like a wave in a homogeneous 

medium. Thus if C + (resp. C ) denotes a cone in R^(d) (resp. R^(d)) 
then Corollary 3.3 and the results of [9] imply 


(5.11) 


^(Ur.C 1 ) - c 2 (±“) p(±«) f ( |p | 2 + q 2 ) |h_(p,q) | 2 dpdq 

Jc* 

■ P(±°°) + |f_(u ± (p,q) ,p,q) | 2 dpdq. 

J r- 


It follows that if 

(5.12) S - {(x,y) : d x < y < d 2 } 
is a slab then 

(5.13) E°°(u f ,S) - 0. 

The Guided Component . Consider the family of cones defined by 

(5.14) C(£,d) - {(x,y) : |y — d| < e|x|} 
where e > 0. Then, in contrast to (5.11), one has 

(5.15) E°°(u^.CCe.d)) - E^.R 3 ) - |A 1 ' 2 1^1J 

for every e > 0 and k ■ 1, 2, 3,••• (see [10, Theorem 5.5]). Finally, 
if S is the slab defined by (5.12) then one can show by means of Theorem 
4.4 that 


*-( V s,-J |f k «\(|pl)’P)l ! 


| ^(y.P) c" 2 (y) p” 1 (y)dy 


dp. 


(5.16) 
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§6. Seal-Infinite and Finite Layers . 

The preceding analysis is extended In this section to the cases 
of semi-infinite and finite layers of stratified fluid. The extensions 
are based on the normal mode expansions for these cases that were derived 
in [12, §9]. Only the principal concepts and results are formulated 
here since the proofs are entirely analogous to those of the preceding 
sections. 

Semi-Infinite Layers . As in [12, §9] the fluid is assumed to 
occupy the domain and to satisfy the Dirichlet or Neumann boundary 
condition. Here the functions p(y) and c(y) are assumed to be 
Lebesgue measurable and satisfy 

(6.1) 0 < p m < p(y) < P M < °°, 0 < c m < c(y) < < 00 

and 

(6.2) |p(y) ~ P(°°)| < C y" a ,|c(y) - c(“)| < C y - ® 

for all y > 0 where P n , p^, p(°°), c m , c^, c(°°), C and a are constants and 

(6.3) a > 3/2. 

As in [12, §9] the acoustic propagators for p(y) and c(y) corresponding 
to the Dirichlet and Neumann conditions will be denoted by A° and A 1 , 
respectively. They are selfadjoint non-negative linear operators in 
*+ - MR£,c _2 (y) p” 1 (y)dxdy) . 

The normal mode functions i^(x,y,p,X) for A^, as defined in 
[12, §9], are parameterized by (p,A) € fl - {(p»A) j X > c 2 (®)jpj 2 }. 
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Their asymptotic form for y ■+■ 00 is given by 


(6.4) 


4,3 ( 


x,y,p,q) 


c(I p 1 J a i(p*x-yq) + R j g i(p*x+yq) 


2tt 


where q - q(|p|,X) - (X c" 2 (°°) - |p| 2 ) 1/5 , c(|p|,X) - (p(®)/4irq( |p| ,X)) 1/2 , 
r3 ■ r3(|p),X) and (R^(|p|»X)| - 1. As in the preceding sections, it 
will be convenient to introduce new parameters: (p,X) -*■ (p,q) 

■ (p»q(|p|»X)) € R^ and normal mode functions 

(6.5) <j>+(x,y,p,q) - (2q) 1/2 c(«) (|/^(s,y,p,X) 

where 


(6.6) X - X(p,q) - c 2 (“)(|p| 2 + q 2 ). 

The asymptotic form of is 


(6.7) <f^.(x,y,p,q) 


c(°°)p ^ (°°) f i(p«x-qy) + R j g i(p*x+qy)\ 

(2ir) Vi \ J 


y ^ oo. 


The second family defined by 


(6.8) <^.(x,y,p,q) - 4>+(x,y,-p,q) 


is also needed. It satisfies 


(6.9) $_(x,y,p,q) ~ 


c(“)p 1/2 («) 

(2ir) W 


j g i(p*x+qy) + R j e i(p«x-qy)^ y 


The expansion theorem of [12, §9] implies that the limits 

M . - 

(6.10) f|(p,q) - L 2 (R?)-lio <fr+(x,y,p,q) f(x,y) c“ 2 (y) p -1 (y)dxdy 

1 1x| <M 







exist. Moreover, if 
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(6.16) v-(t,x,y) - (t^(x,y,p,q) exp (-It u(p,q)) h^(p,q)dpdq 

1 Jn 3 


where w(p,q) * c(°°)>/|p | 2 + q 2 . Moreover, as in §3 one can write 


(6.17) 4*^(x,y,p,q) 


c(“)p l/2 (») 


(2ir) 


3/2 


|.1(P"«W) I^(y,p,q) +e 1(p ** -qy) R 


where 


(6.18) 


lim 1^(y,p,q) - 1, 
Y** 


lim R^(y,p,q) - R^(p,X). 

Y** 

Then, proceeding as in §3, one can prove the following analogue of 
Theorem 3.1. 

Theorem 6.1 . For every h 6 3f + let v£(t,’) be defined by 
(6.19) Vj(t.x.y) - exp (-itc(“) Aj /2 ) h„(x,y), (x,y) € Rj. 

where h 0 € L 2 (R 3 ) is the function whose Fourier transform is 


( 6 . 20 ) 


Then 


b<j(p»q) - 4 


c (oo) p ^2 (») hJ(p,q), (p,q) e R®, 

0 , (p»q) € R*. 


(6.21) 


lim |v f (t,«) - v|(t,»)|„ - 0. 


i!(y>p,q) 
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Transient Guided Waves . For both semi-infinite and finite 
layers the form of the guided components v^(t,•) is precisely the same 
as for the case of an infinite layer. Thus the analysis of §4 applies 
unchanged to these cases. 
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